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The nuclear stopping, the elliptic flow, and the HBT interferometry are calculated by the UrQMD 
transport model, in which potentials for "pre-formed" particles (string fragments) from color flux- 
tube fragmentation as well as for confined particles are considered. This description provides stronger 
pressure at the early stage and describes these observables better than the default cascade mode 
(where the "pre-formed" particles from string fragmentation are treated to be free-streaming). It 
should be stressed that the inclusion of potential interactions pushes down the calculated HBT 
radius Ro and pulls up the Rs so that the HBT time-related puzzle disappears throughout the 
energies from AGS, SPS, to RHIC. 
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One of the most important discoveries from SPS and 
RHIC experiments in the beginning of this century is that 
a new form of matter is produced in high energy ion colli- 
sions that seems to behave as a strongly interacting quark 
gluon plasma (sQGP) jll,|l,[3,H. Furthermore, at RHIC 
energies, the sQGP seems to be a nearly "perfect liquid" 
with small viscosity [1, Q. Although hydrodynamical 
studies are a useful tool for the theoretical investigations 
at RHIC energies, hydrodynamics fails at lower energies. 
Therefore, a relativistic dynamic transport approach, in 
which an entire equilibrium is not pre-assumed, is nec- 
essary to explore the excitation functions of various ob- 
servables. In addition, the transport model has the ad- 
vantage of observing the whole phase-space evolution of 
all particles involved in a microscopic fashion. 

In order to detect this new matter - quark gluon 
plasma (QGP) - many theoretical suggested observables 
have been discussed and argued frequently. Among 
these, the Hanbury-Brown-Twiss interferometry (HBT) 
or Femtoscopy technique has been used widely to extract 
the spatio-temporal information of the particle freeze-out 
source. It would be a vital discovery if there is an nontriv- 
ial transition in the spatio-temporal characteristics of the 
source, when going from low to high beam energies 
Experimentally, the HBT parameters have been scanned 
thoroughly over the energies from SIS, AGS, SPS, up to 
RHIC, unfortunately, no obvious discontinuities do ap- 
pear [8]. 

However, this "null" result does not mean that the 
HBT technique comes to an end. Non-Gaussian effects 
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might shadow the possible energy dependence of the HBT 
parameters E-puzzle" ) , and this to pic has been quickly 
improving in the recent years BEini- Secondly, even 
with a Gaussian parametrization, we found that the HBT 
time-related puzzle (^'t-puzzle'^) is present at almost all 
energies from AGS to RHIC [I4I . In order to understand 
the origins of these "puzzles" , it is necessary to dig deeper 
into the dynamics of the heavy ion collisions. 

The hydrodynamic model, in which various kinds of 
equations of state (EoS) with latent heats are considered, 
successfully explained the elliptic flow V2 at transverse 
momentum pt < 2GeV/c at RHIC energies [H [3 [H. 
However, hydrodynamics failed to explain the experimen- 
tal Ro/Rs data {f§\ and the references therein), i.e. the 
ratio between the HBT radii in outward and sideward 
directions, where the "t-puzzle" initially occurred. But 
recently, it has been realized that a direct (fully "apple- 
to-apple" ) comparison between hydro calculations and 
data is required 16]. 

In the transport models, the effective EoS is softened 
by the existence of resonances, but the re-scattering be- 
tween particles will reinforce to some extent the stiff- 
ness of the dense matter. At SPS and RHIC energies, 
the color-string excitation and fragmentation is further 
considered for particle production and subsequent evolu- 
tion in the ultra- relativistic/relativistic quantum molecu- 
lar dynamics (UrQMD /RQMD) hadron-string transport 
models. The formation time of the hadrons from string 
fragmentation is determined by a "yo-yo" mode [17, 1§|. 
During the formation time, the "pre-formed" particles 
(string fragments that will be projected onto hadron 
states later on) are usually treated to be free-streaming, 
while reduced cross sections are only included for leading 
hadrons. The idea of "pre-formed" hadrons is a well- 
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known concept in the description of deep inelastic scat- 
tering data (DIS) [l^ll^]- It has also been found that the 
"pre-formed" hadrons from the color fluxtube fragmen- 
tation behave already very hadronlike, before the hadron 
formation is fully completed [21,]. In previous calcula- 
tions [13, [H, [HI, the interaction of the newly created 
"pre-formed" particles from the string fragmentation is 
usually not taken into account for simplicity. The facts 
have shown that the effective EoS at the early stage in 
the past model simulations is too soft and the early pres- 
sure is inadequate to match the calculated elliptic flow 
to the experimentally observed flows at RHIC energies 
[23I [23 | . In order to make up for this discrepancy in the 
early-stage EoS, it is essential to modify the dynamics 
of the early stage with/ without new degrees of freedom, 
such as done in a multi-phase tr ansp ort (AMPT) [ll| . 
the hadronic rescattering (HRM) [2^, the quark molec- 
ular dynamics (qMD) [2y| , and the parton hadron string 
dynamics (PHSD) ^ model. 

Besides the new matter which is dominant at the early 
stage, the potentials between confined hadrons which 
play a role during almost all time of the collision process 
deserve attention as well. In the previously low energy 
UrQMD calculations, the Hamiltonian H served for the 
EoS consists of the kinetic energy T and the effective 
two-body interaction potential energy U. The poten- 
tial energy U is indispensable, which includes the two- 
body and three-body (which can be approximately writ- 
ten in the form of two-body interaction) Skyrme- (also 
called as the densit y d ependent terms), Yukawa-, and 
Coulomb-terms [13, 13- Recently, in order to be more 
successfully applied in the intermediate energy region 
{El, < 2A GeV) , more terms are considered ^2§„ 23} , those 
are, the density-dependent symmetry potential (essen- 
tial for isospin-asymmetric reactions at intermediate and 
low energies) and the momentum-dependent term. At 
higher energies, i.e., AGS/FAIR and SPS, the Yukawa- 
and symmetry- potentials of baryons becomes negligi- 
ble, while the Skyrme- and the momentum-dependent 
part of potentials still influence the whole dynamic pro- 
cess of HICs though weakly. In Ref. [lo], the effects of 
the mean field with momentum dependence on collective 
flows from HICs at 2-158A GeV energies were studied 
by a Jet AA Microscopic Transportation Model (JAM) 
and it has been found that the momentum dependence in 
the nuclear mean field is important for the understand- 
ing of the proton collective flows at AGS and even at SPS 
energies. Further, by adopting the same soft EoS with 
momentum dependence (SM-EoS) on confined baryons, 
we also found that the calculated Rq / Rs ratio from HICs 
at AGS energies can be vastly improved when compared 
to the measured data [H, jsif. 

In this work, we will follow this idea and investigate the 
HBT radii at higher beam energies, i.e. SPS and RHIC. 
Furthermore, the other observables, such as the rapidity 
distributions of net protons (p — p) and pions, and the el- 
liptic flows of charged hadrons, are checked at the same 
time. As an initial attempt, we consider potentials for 



both formed and "pre-formed" particles: 1), for formed 
baryons, the SM-EoS [13] is chosen. The Coulomb poten- 
tial is also considered, as in the low-energy calculations 
[H, dH- 2), for formed mesons, only the Coulomb po- 
tential is optional in calculations, while other potentials 
are not considered as in previous calculations at low en- 
ergies [2^ [m • 3) , for "pre-formed" particles from string 
fragmentation, the similar density dependent (Skyrme- 
like) terms [43] as the formed baryons are used, but 
without the Yukawa, the Coulomb, and the momentum 
dependent terms. 4), the "pre-formed" mesons act like 
"pre-formed" baryons but with a reduction factor (2/3) 
due to the quark-number difference. Furthermore, 5), 
there is no potential interaction between "pre-formed" 
baryons and formed baryons so far. Correspondingly, 
the hadronic density is calculated by ph = J2j^i^i^jPij^ 
where Cij ~ 1 for Baryons, 2/3 for "pre-formed" mesons, 
and for formed mesons. Pp is a Gaussian in coordi- 
nate space. As stated in [sO, HH, the relativistic effect 
(Lorentz transformation) on the relative distance and rel- 
ative momentum between two particles i and j has been 
considered. 

To set the stage, Fig. [T] shows the rapidity dependence 
of the multiplicities of net protons (top) and negatively 
charged pions (bottom) for central (< 5% of total cross 
section or) Pb-|-Pb coUisions at Eb = 158A GeV. The 
experimental data (solid stars) are taken from [3^ . [33j , 
and the open stars represent the data reflected at mid- 
rapidity. The dash-dot-dotted lines are the cascade mode 
calculations, the dashed lines represent the calculations 
with SM-EoS of confined baryons, while the solid lines il- 
lustrate the results with the additional density-dependent 
potentials for "pre-formed" particles. While the cas- 
cade calculations give a reasonable rapidity distribution 
for the pions, the rapidity distribution of net protons 
is Gaussian - as seen in previous UrQMD calculations 
[2^ . and not fully consistent with the data (which has a 
two-bump structure at ?/ ~ ±1.2). After considering the 
baryonic potentials, the distribution only changes slightly 
at projectile and target rapidities. However, when the 
density-dependent EoS for "pre-formed" particles is addi- 
tionally taken into account, one observes that two bumps 
in the rapidity distribution of net protons appear. This 
implies that, due to the interactions of "pre-formed" par- 
ticles at early stage, the additional pressure increases the 
effective stopping power as well as the central density at 
the compression stage. The total pion multiplicity re- 
mains basically unchanged on a level below 5% at this 
beam energy. Thus, the current pion distribution with 
potentials is acceptable for further investigations. 

Using the analyzing program "correlation after- 
burner" (CRAB y3.0/3), which is contributed by S. Pratt 
[33 , [35} . we calculate the HBT correlator and further 
fit it as a three-dimensional Gaussian form under the 
Pratt convention, i.e., the longitudinally comoving sys- 
tem (LCMS), which is expressed as 
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FIG. 1: Rapidity distributions of net protons (top) and neg- 
atively charged pions (bottom) from central (< 5% of total 
cross section) Pb+Pb collisions at Ei, — 158A GeV. The ex- 
perimental data (solid stars) are taken from [s^, H^l, while 
the open stars represent the data reflected at j/ = in the 
center-of-mass system. The dash-dot-dotted lines are the cal- 
culations with cascade mode, the dashed lines represent the 
calculations with potentials for confined baryons (SM-EoS), 
while in the solid lines, the density-dependent potentials for 
"pre- formed" particles are further taken into account. 



C{qo,qs,qL) = l+Xe^p{-Rlql-Rl,ql,~Rlql-2Rl,^qoqL) 

(1) 

In Eq. ^ the A is the incoherence factor, Rl, Rq, and 
Rs are the Pratt radii in longitudinal, outward, and 
sideward directions, while the cross-term Rql pla-ys a 
role at large rapidities, qi is the pair relative momen- 
tum in the i direction. Fig. [2] shows the transverse mo- 
mentum kx dependence (ky = (Pit + P2t)/2) of the 
HBT-radii Rl, Roi *md Rs (at midrapidity) for central 
HICs at AGS {Et = 2A, 4A, 6A, and 8A GeV), SPS 
{Eb = 20 A, 30A, 40A, 80A, and 158A GeV), and RHIC 
i^/sNN = 130 and 200 GeV) energies. The physical cuts 
used are same as those listed in |12|. The data are from 
the E895^^NA49j,^STAR, and PHENIX collaborations 
dE [13, iE m, m. The calculations with cascade 
mode are shown by lines with down-triangles while the 
calculations with various potential treatments are shown 
with other lines with different symbols. At Et, = 2A GeV, 
we show a comparison of calculations between with and 
without Coulomb potentials for confined mesons. It is 
seen that the two-body mesonic Coulomb potential be- 



fore freeze-out affects the HBT radii only very weakly at 
about kx < lOOMeV/c, and the ratio between Rq and 
Rs values is not altered consequently. At AGS energies, 
besides the nuclear potentials, the Coulomb potential of 
confined mesons is also considered (solid lines with cir- 
cles), while at SPS and RHIC energies, the Coulomb po- 
tential for confined mesons is switched off in order to 
save computing times and to avoid problems with non- 
local interactions. It is known that, below AGS energies 
hadrons are dominantly produced from the decay of reso- 
nances so that the considered potentials for "pre-formed" 
particles from string fragmentations have no effect on the 
HBT radii [13, [13 • At higher energies "pre-formed" par- 
ticles play a dominant role for the early stage dynamics, 
which can be seen, e.g., by the comparison of solid lines 
and the lines with triangles at y^SNN — 200 GeV. So far 
the question of a steeper fcT-dependence of HBT radii 
than data at high SPS and RHIC energies, similar to the 
AMPT transport model calculations with parton scatter- 
ings is still open. 

The most exciting results show up in the transverse 
space. Roughly speaking, the nuclear potentials lead to 
a smaller Rq (especially at large fcr) but a larger Rs 
(especially at small kx) generally allowing for a better 
description of the kx dependent data. Certainly, our 
present crude approach is only semi-realistic as a descrip- 
tion of the dynamic processes at the early stage of the 
reaction. However, one can still come to the conclusion 
that the potentials help to obtain a better description of 
the pion freeze-out in the transverse space. 

In Fig. [3] we depict the excitation function of the 
Ro/Rs ratio at small kx'- at AGS and SPS energies, the 
results at kx 150 ± 50MeV/c are shown, while at RHIC 
energies, the kx is set to 200 ± 50MeV/c. The exper- 
imental data within the transverse momentum regions 
are compared with the calculations with and without po- 
tentials. As seen before [H, [11, [1E| , in the cascade mode 
.(dash-dot-dotted lines with open rectangles) the Rq/Rs 
ratio, which is related to the duration of the emission t 
(t ~ \/ Rq ^ R^)i is larger than the experimentally ob- 
served values at all investigated energies. Furthermore, 
the increase of the Rq /Rs ratio with beam energies at 
AGS is not seen in these calculations. When the SM- 
EoS is considered for confined baryons (dashed lines with 
open up-triangles), the Rq/Rs ratio decreases compared 
to the result in the cascade mode and reproduces the en- 
ergy dependence of the data up to the lower SPS energies. 
At high SPS and RHIC energies, however, the Rq/Rs ra- 
tio increases further with increasing beam energies and 
deviates strongly from the data at top RHIC energies. 
Here, the ratio nearly approaches the value obtained in 
the cascade mode. It implies that the potential of con- 
fined hadrons is increasingly losing its importance with 
increasing beam energy. At SPS and RHIC energies the 
deviation from data can be interpreted by the absence of 
interactions of "pre-formed" particles from string frag- 
mentation since at Ef, ^ 30 — 40A GeV the measured 
Ro/Rs ratio begins to flatten/decrease in contrast to 
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FIG. 2: Transverse momentum kr dependence of the HBT radii Rl, Ro, and Rs (at midrapidity) for central HICs at AGS 
{Eb = 2A, 4A, 6A, and 8A GeV), SPS {Eb = 20A, 30A, 40A, 80A, and 158A GeV), and RHIC {^fs^N = 130 and 200 
GeV) energies. The data are indicated by sohd stars, which are from E895, NA49, STAR, and PHENIX collaborations 
[3^. IstI. [3^. Isgl. [iol. I4H ■ The calculations with cascade mode are shown by lines with down-triangles. At AGS energies, the 
calculations with potentials for confined baryons (SM-EoS) as well as a Coulomb potential for confined mesons are shown by 
the lines with circles. At AGS E^ = 2A GeV, SPS, and RHIC energies, the results with potentials for both "pre-formed" 
particles and confined baryons but without Coulomb potential for mesons are shown by lines, while the lines with up-triangles 
correspond to the calculations without any potentials for "pre-formed" particles but with potentials for confined baryons. 



the aforementioned calculations. However, with the con- 
sideration of a density dependent EoS for "pre-formed" 
hadrons (solid lines with open down-triangles) this be- 
havior can be reproduced in the calculations. With this 
approach, at RHIC energies, the ratio Rq/Rs is about 
unity (slightly below data) which implies that the com- 
pressibility provided by the currently employed potential 
is too strong. Nevertheless, since we are not aiming at 
fitting the whole excitation function of the data, we leave 
this detail aside for further detailed explorations. 

Therefore, our explanation of the HBT time related 
puzzle is different from the investigations in Ref. (46j 
where Cramer, et al., considered a single-pion optical 
potential which simulates its interactions with the dense 
medium. While in this work, the microscopic "nuclear" 
potentials are considered for "pre-formed" hadrons from 
the string fragmentation. The Hamilton's equations of 



motion of the particles represented by Gaussian wave 
packets are solved microscopically. After that, the freeze- 
out phase space is analyzed numerically by the CRAB 
program, in which the plain and distorted (by Coulomb 
and nuclear potentials after freeze-out) wave functions 
are optional. 

Let us finally check the modification of the early trans- 
verse expansion by a look into the elliptic flow. In the up- 
per plots of Fig. [4] we present the transverse momentum 
Pt dependence of the elhptic flow V2 {v2 — {{Px ~'Py)/Pt)j 
where p^ = p"^ + p^) of charged pions from Au-|-Au col- 
lisions at RHIC y/SNN = 200 GeV. Two sets of cen- 
tralities are chosen: 5% < a/(TT < 10% (in (a)) and 
20% < ct/ctt < 30% ((b)). To avoid non-flow effects [13, 
the four-particle integral cumulant V2 data (f2{4}) are 
shown with stars. Similar to those in Fig. [31 the UrQMD 
calculations without and with potentials are illustrated 
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FIG. 3: Excitation function of the Ro/Rs ratio at small kr- 
The data are indicated by solid stars. The dash-dot-dotted 
lines with open rectangles are results under cascade mode, the 
dashed lines with open up-triangles represent the calculations 
with potentials for confined baryons, while the solid lines with 
open down-triangles are results with potentials for both "pre- 
formed" and confined particles. 
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for comparison with data. It is clear that due to the ab- 
sence of early pressure, the elliptic flow in the cascade 
mode (dash-dot-dotted lines with open rectangles) can 
not grow as fast as the data with the increase of trans- 
verse momentum in both centralities. After switching on 
the SM-EoS for confined baryons (dashed lines with open 
up-triangles), the flow is seen to increase slightly espe- 
cially at high transverse momenta. With the density de- 
pendent potentials for "pre-formed" particles (solid lines 
with open down-triangles) , the elliptic flow at large trans- 
verse momenta rises further, however, at small transverse 
momenta it is driven down and deviates from data espe- 
cially at large centralities. This phenomenon calls for a 
more complete description of the dynamics of HICs at 
high energies in which the (in)elastic scatterings of par- 
ticles (confined hadrons and deconfined partons) should 
be considered in more detail. To make a whole com- 
parison, we show in Fig. 2] (c) the centrality depen- 
dence of the ratio between different theoretical V2 values 
of charged hadrons and the four-particle cumulant data 
{"v2^ / V2^^ {4:}" , here the flows are integrated over pt and 
rj). As stated before, the flow obtained from the cas- 
cade calculation is only about 55% of the experimental 
one. The contribution of potentials of confined baryons 
to flow is only at the level of about 4%. While the poten- 
tials of "pre-formed" hadrons increase the flow in more 
central collisions. At large centrality, the potential effects 
are relatively weak. 

In summary, besides a soft equation of state with mo- 
mentum dependence, which is required especially at low 
energies, a density dependent (Skyrme-like) potential has 
been tested for the "pre-formed" particles from string 
fragmentation in the UrQMD transport model. Although 
the form of the potential for the new phase is simple and 



FIG. 4: Top plots ((a) and (b)): Transverse momentum pt de- 
pendence of the elliptic fiow V2 for Au+Au collisions at RHIC 
y^sjvjv = 200 GeV. The pseudo-rapidity cut is I77I < 1.3. Two 
sets of centralities are chosen: 5% < cr/ar < 10% (in (a)) and 
20% < (t/(tt < 30% ((b)). Stars represent the four-particle 
integral cumulant 112 data (w2{4}). Dash-dot-dotted lines with 
open rectangles represent the cascade calculations. Dashed 
lines with open up-triangles are the results with confined bary- 
onic potentials, while the solid lines with open down-triangles 
are the results with potentials for both "pre-formed" particles 
and confined baryons. Bottom plot ((c)): Centrality depen- 
dence of the ratio between theoretical and experimental fiows 
of charged hadrons, which are integrated over pt and rj. 

rough (ideally, the EoS of the new phase should be based 
on the first-principle lattice QCD calculations [3), it 
provides new insights into the dynamics of HICs. It re- 
produces the proper stopping power for net protons. It 
also improves the elliptic flow at large transverse mo- 
menta. Most importantly however, it decreases the HBT 
radius Rq but increases the Rs so that the experimental 
duration time related quantity — here the Rq / Rs ratio 
— can be reasonably reproduced throughout the ener- 
gies from AGS, SPS, up to RHIC. Although some open 
problems remain due to the simple treatment of the in- 
teraction of the particles at the early stage. 
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